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ABSTRACT 



The principle of virtual work v/as used to determine 
several typical difference equations governing the dis- 
placements of a taut string subjected to a distributed 
transverse load. The integrals for the strain energy and 
the external work were approximated by both the trapezoidal 
rule and Simpson's rule. Further^ the integrand in the 
strain energy expression was replaced by two different 
finite-difference approximations. For comparison;, two 
typical difference equations v;ere developed from the 
differential equation for the string using central finite 
.a-i -Ffp-rprsopc; and the Hermitian differencing scheme. The 
finite difference equations j all of i^hich are tridiagonalj 
were evaluated by comparing truncation error. 
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I . INTRODUCTION 



The principle of virtual work is a very useful tool 
for finding the displacements of an elastic body produced 
by given forces In this thesis the principle is 

used to determine several finite difference equations 
governing the displacements of a taut elastic string sub- 
jected to a distributed transverse load per unit lengthy 
q. The strain energy of the string, U, is given by 




in X'/hich F is the constant internal force in the stringy 
w is the displacement L is the length, and x is the axial 
coordinate. The positive directions of F, w, q, ctud 
are indicated in Figure 1. 




Figure 1. A STRING SUBJECTED TO A TRANVERSE LOAD 

AND A TENSILE FORCE 
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The principle of virtual v;ork is 



<^U 



( 1 . 2 ) 



where is the variation of U with respect to w and /w 
is the work performed by q moving through the virtual 
displacement <Jw. Hence 




(1-3) 



The trapezoidal rule and Simpson's rule are used to 
numerically evaluate the continuous integrals in Equations 
1.1 and 1.3 when the string is divided into k intervals of 
length h. The displacements of the string are defined at 
the end points or nodes of these intervals ^ and the 
continuous first derivative in U is replaced by a finite 
difference expression. Several approxiimations to the first 
derviative are considered. This procedure produces a 
lengthy equation in terms of tne aispiacemen is and loads 
at the finite number of nodes . Taking the variation of this 
equation with respect to the displacement at each of the 
nodes^j in accordance with the principle of virtual work, 
leads to a set of simultaneous linear algebraic equations. 
There is one such set of equations for each numerical scheme. 
Within each set there are one or two typical difference 
equations. The schemes are evaluated by examining the 
truncation errors associated with each typical equation. 

^The variation must satisfy the geometric boundary 
conditions, i.e.,^w=oatx=o and x = L. 
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A second procedure for finding the displacements of 
an elastic body is to use a finite difference approximation 
of the governing differential equation. For the string, 
the differential equation is 



Substituting a finite difference approximation for the 
second derivative leads to a set of simultaneous linear 
algebraic equations of tridiagonal form. Tv;o typical 
finite difference equations are derived and compared with 
the typical equations developed from the virtual vfork 
principle . 



F 




+ q = 0 



( 1 . 4 ) 
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II. DEVELOPMENT OF THE FINITE DIFFERENCE EQUATIONS 



USING THE PRINCIPLE OF VIRTUAL WORK 



A. THE TRAPEZOIDAL RULE A^TD SIMPSON'S RULE 

There are many formulas available for numerically 
evaluating integrals. The trapezoidal rule and Simpson's 
rulej two of the most commonly used^ v/ere chosen for this 
thesis. Both of these find approximations to integrals 
by passing an interpolating function through several nodes 
and integrating the interpolating function. The trape- 
zoidal rule uses a straight line as the interpolating 
function between x and x , ^ , where x„ is the value of 
X at the th node_, and is given as ['2] 




± } 



s+1 



hr '1 h" 

f(x)dx= 2|^f(Xg_^^) - - j2 



r dLin 



dx 



2 ” 



in which x < T<x The last term in Equation 2.1 

s _ L. - s+1 

is knovm as the truncation error. The integral from x=0 
to x=L is 




( 2 . 2 ) 

f(x)dx = hjl-f (0)+f (x^)+. . .+f (x^_^)+^f (L)J 



- h^L 

12 dx-^ 



in which 0 ~ TT — L 

In Simpson's rule a parabola is used as the interpolat- 
ing fxmction passing through the nodes x , ^ ^ x . 

s — 1 s s +x 
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The integral from x=^0 to x=L is 

CL 

f(x)dx = 

0 3 



(2.3) 






f(x)dx = h j^f(0)-i4f(x3^)+2f(x2) + . . .+2f{x^_^) 

] - A dVm 



90 dx 

where 0 5 ZT - L and k is even. 



B. DIFFERENCE APPROXIMATIONS TO THE FIRST DERIVATIVE 



The equation containing the conventional central 
finite difference approximation to the first derivative 
at Xg and the truncation error is [2] 






w _ ) 
s-1 ^ 






(2.4) 



in which x ^ < Vs ^ x 

s-1- s+1 



Consider the scheme in which the displacements are 
defined at the nodes and the first derivatives are defined 
at the center of the intervals as shown in Figure 2. This 
scheme is referred to as the staggered differencing scheme 
or half-station scheme. The central finite difference 
equation for the staggered scheme is derived in Appendix 
and is 




= 1 (w 
s+i ^ 



s+1 



- w^ ) - h^ (d^w) 

= TO ^ 



(S-5) 



where x T , i 5. x ,t. Hereafter, unless noted otherv;ise 
s ^s+2 s+1 ^ 

Xn denotes that X is located in the interval x ^ to x 
s s-1 s+1 

and2^g+i is in the interval x^ to and T is in the 

interval x=0 to x=L. 
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Figure 2. THE STAGGERED SCHEME 



Another finite difference equation used to replace the 
first derivative at x is this higher order expression 

o 



/ _dw\ -w )+Il 

I dx/ ^ 2h ^ s+1 s-1^ 12 F 



(q. t) + /d^w\ 

^^s+1 ^s-1^ 

( 2 . 6 ) 

The corresponding higher order staggered scheme equation is 

/dw) = 1 (w (q ^ - Q ) q 71^^ / d^wl 

h s+1 si 2ifp ' s+1 *13^ 5^ 



(2-7) 



-4-i 



Equations 2.6 and 2.7 are derived in Appendix A. 



C. DEVELOPMENT OF THE TYPICAL FINITE DIFFERENCE EQUATIONS 
The first typical difference equation is derived using 
the trapezoidal rule and the central finite difference ex- 
pression. Applying Equations 2.2 and 2.4 to Equation 1.1 
leads to: 
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2 

h FL 
24 



£h ('‘s+2- «s ) 



dx 



dw\^ 

'^Jz 

2 



( 2 . 8 ) 



Only those terms in which w appears are listed in Equation 

o 

2.8 because only the typical equation is of interest. 

Taking the variation of Equation 2.8 with respect to the 



displacement leads to 



^ - P (w3_2-2w +w )-hiP r 

iWs 4h ^ 12 I 



(2-9) 

/d\rj 



s+1 



Taking the variation of the external work with respect to 

w leads to 
s 



= hq 



( 2 . 10 ) 



^ This equation does not take into 

s-2 s+2 ^ 



Note in Equation 2.9 that the displacements are required 

at ^ and 
s-2 

consideration the fact that the displacements at and 

are available. Hence ^ it appears reasonable to use 

these displacements instead of those at and 

s-2 s+2 

Consequently^ let the interval h be reduced by a factor of 
one half and change the subscripts accordingly to give a 
set of equations that do not omit any nodes. The modified 
equations are : 
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( 




( 2 . 11 ) 



iv _ h 

^''s 2 



( 2 . 12 ) 



It can be shown that equations 2.11 and 2.12 are identical 
to the equations developed using the trapezoidal rule and 
the staggered difference scheme. Using the staggered 
scheme eliminates the necessity to halve the internal. 

As an example of the derivation of the higher order 
equations, a typical equation is derived using Simpson's 
rule in conJunr*t"' on i*n‘ th the higher order staggered scheme. 
The equation for U based on Simpson's rule can be given in 
the form 




2 





J 



4 

h L d 



90 



(2.13a) 



Equation 2.13a could also have been given as 






(2.13b) 
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The typical equations associated with Equations 2.13a and 
2.13h are different. The equation derived from Equation 
2.13a is referred to as the 2-4-2 scheme and the equation 
derived from Equation 2.13h is referred to as the 4-2-4 
scheme . 

Substituting Equation 2.7 into Equation 2.13a gives 




Taking the variation of U with respect to w leads to 



<Jw 



2F 

3h 



+ 2«s+l 



) “ (<1 - 3q + 2q ) 

^ 36 ^ s-1 s s+l' 



+ 



7h^F 

8640 




(2-15) 



The variation of the external work with respect to w^ for 
the integral evaluated by Simpson 's rule with the 2-4-2 
scheme yields 



<£w _ 4h 

Wg 3 



(2.16) 
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All of the typical difference equations that were derived 
using the principle of virtual vrork are given in Tables I 
and II. Table I contains the equations derived from the 
trapezoidal rule and Table II contains the equations derived 
from Simpson's rule. 

The trimcation errors listed in Tables I and II are only 
part of the total truncation error. These trimcation errors 
are associated with the finite difference evaluation of the 
derivatives in the strain energy integral. Note that they 
consist of the difference of two derivatives. As a con- 
sequence, they can be converted to a higher order of h. As 
an example, from the first equation of Table I 
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TRAPEZOIDAL RULE FINITE D DEFERENCE EQUATIONS 
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N - non-staggered interval; S - staggered interval; 2-4-2 - Simpson's rule^ 2-4-2 scheme; 
4-2-4 - Simnson's rule, 4-2-4 scheme 
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III. THE DIFFERENTIAL EQUATION APPROACH 



A. CENTRAL FINITE DIFFERENCE FORMULATION 



The differential equation for the string is given by 
Equation 1.4. Using the conventional central finite difference 
expression to evaluate the second derivative at s leads to 



C 2] 

F(w -2w +w )/h 
^ s-1 s s+1^ 




(3.1) 



B. HERMITIAN FINITE DIFFERENCE FORMULATION 

The Hermitian formula for the second derivative^ as 
derived in Appendix Cj is 




Applying Equation 1.4 at 
gives 



s 



and to Equation 3.2 

s+1 



F (w - 2w + w ) + 1 f q 



s+l‘ 



12 



^s-1 



+ lOq 



+ q ) 



+ !LZ ( 

240 I dj^/Tg 



0 



(3.3) 



Equation 3.3 is the typical difference equation for the 
Hermitian finite difference scheme. 
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Equations 3.1 and 3.3;. the t\^o typical difference 
equations, are listed in Table III. 
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FINITE DIFFERENCE APPROXIMATION TO TEE STRING'S DIFFERENTIAL EQUATION 




22 



IV. CONCLUSIONS 



All of the equations derived give tridiagonal matrices 
centered about the principal diagonal. For a given number 
of nodes a computer program using these equations should 
have approximately equal storage requirements computation 
time and round-off error. 

There are tvro ways to decrease the truncation error. 

One is to use higher order difference equations and the 
other is to reduce the interval^ making more nodes. Normally, 
both of these cause greater round-off error, storage 
requirements, and computation time, since the first increases 
the number of non-zero terms in the matrix and the second 
increases the size of the matrix. 

The higher order equations derived for this thesis do 
not increase the size of the displacement matrix but do 
change the load matrix from a diagonal to a tridiagonal 
matrix. 

The object of any numerical analysis method is to obtain 
as accurate an ansvrer as possible within the shortest com- 
putation time. Of the two finite difference solutions to 
the differential equation, the Hermitian method is the better 
of the two from the standpoint of truncation error. 

For the equations derived from the energy solution, 
the truncation error is composed of three parts. Two parts 
are associated with the evaluation of the strain energy and 
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work integrals and are of O(h^) for the trapezoidal rule 

4 

and 0(h ) for Simpson's rule. The third is the one 

associated v;ith the evaluation of the derivatives in the 

strain energy integral and is listed in Tables I and II. 

2 

This part v;as of 0(h ) x^rhen the central finite difference 

4 

expression was used and of 0(h ) when the higher order 
finite difference expression was used. 

From the comparison of the equations made in this 
thesis by partial truncation error it would appear that the 
energy equations derived using Simpson's rule and the 
higher-order finite difference equations are the better 
equations. However_, it must be stressed that this evaluation 
is based only on the apparent trutication error and is, 
therefore, not conclusive. 

There is some difficulty in establishing the total 
truncation error associated with the equations derived from 
the energy approach. Consequently, it is recommended that 
these equations be evaluated on the computer so that a com- 
plete comparison can be made not only of truncation error 
but also of storage requirements and computation time. 
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APPENDIX A 



THE HIGHER ORDER FINITE DIFFERENCE EQUATION FOR THE FIRST 

DERIVATIVE 



V/riting the Taylor series for taking the 

derivative twice gives 




(A.l) 

and for w , gives 
S "" _L 




(A. 2) 

Subtracting Equation A. 2 from Equation A.l gives 




(A. 4) 
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Substituting Equation A ,4 into the central finite 
difference equation. Equation 2.4, for the non-staggered 
interval gives 




2h 



' s+1 s-1^ 





(A. 5) 



The differential equation for the string is 




(A. 6) 



Substituting Equation A. 6 into Equation A. 5 and truncating 
gives 




(w. 



I n "'W 

S+1 S 



-4 



s+1 %-l 



— 4 / .h \ 
) + III I'il' 



36oldx5 






(A. 7) 



Equation A. 7 is the higher order finite difference equation 
for the non-staggered scheme. 
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APPENDIX B 



THE DERIVATION OF THE FORMULAS AND EQUATIONS USED WITH 
THE STAGGERED SCHEME 

A. THE CENTRAL FINITE DIFFERENCE EQUATION FOR THE FIRST 
DERIVATIVE 

The Taylor series for w is 

s+l 




Subtracting Equation B.2 from Equation B.l gives 




(B.3) 
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gives 




(BA) 

By truncating after the term in w the central finite 
difference equation for the first derivative for the 
staggered interval is 




where x T. i - x . (B.5) 

s '^s+t s+1 ^ 

B. THE HIGHER ORDER ETNTTE DTEFF.RRNCE EO.UATION FOR FIRST 
DERIA^ATIVE 

Writing the Taylor series for v/ and taking the 

s+1 

derivative twice gives 




and for w have 
s 
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1 











(B.7) 



Subtracting Equation B.7 from Equation B.6 gives 




(B.9) 



Substituting Equation Be 9 into Equation B.4 gives 




Truncating the series gives 
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The differential equation for the string is 




Substituting this into Equation B. 11 gives 




Equation B.12 is the higher order finite difference 
equation for the staggered scheme. 
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APPENDIX C 



THE DERIVATION OF THE HERMITIAN FINITE DIFFERENCE EQUATION 

FOR THE SECOITO DERIVATIVE 
Adding Equations A.l and A . 2 gives 



(S) ;#) ■"‘(‘M 



+ 



iT 

12 I dx 



d) 






( 0 . 1 ) 



Solving Equation C.l for 






/aa iTaan . AdA') , 7dA'\ 1- 1£ 

[dxy s 7 (dx-/ I ax / s ^ dx /3+j *-\ax,g 

(0.2) 

The Taylor series derived expression for the second 
derivative [iJ is 

/df«] = _| (w -2w +w ) - + . . 

(dx^y^ s-l s s+1^ 12 ( dx/ g 360 I /s 



(C.3) 
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1 






Substituting Equation C.2 into Equation C.3 gives 




— o (w„ , -2w +w 
s-1 s 



S+1' 




s-1 



-2 





(C.4) 



Solving Equation C.4 for (w -2w +w ) gives 

S ~ _L S S "HI 




Equation C,5 is the Herrnitian equation for the second 
derivative . 
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